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Abstract 

We present results of near-infrared photometric and spectroscopic observations of 
mid-infrared (MIR) sources that dramatically brightened. Using IRAS, AKARI, and 
WISE point source catalogs, we found that 4 sources (IRAS 19574+491, V2494 Cyg, 
IRAS 22343+7501, and V583 Cas) significantly brightened at MIR wavelengths 
over the 20-30 years of difference in observing times. Little is known about these 
sources except V2494 Cyg, which is considered a FU Orionis star. Our observa¬ 
tion clearly resolves IRAS 22343+7501 into 4 stars (2MASS J22352345+7517076, 
2MASS J22352442+7517037, [RD95] C, and 2MASS J22352497+7517113) and 
first JHK S photometric data for all 4 sources are obtained. Two of these stars 
(2MASS J22352442+7517037 and 2MASS J22352497+7517113) are known as T Tau 
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stars. Our spectroscopic observation reveals that IRAS 19574+9441 is an M-type 
evolved star and V583 Cas is a carbon star. 2MASS J22352345+7517076 is probably 
a YSO, judging from our observation that it has featureless near-infrared (NIR) spec¬ 
trum and also showed dramatic brightening in NIR (about 4 magnitudes in K s -band). 
The possible reasons for dramatic brightening in MIR are discussed in this paper. 

Key words: infrared: stars — (stars:) circumstellar matter — stars: late-type 
— stars: pre-main sequence 


1. Introduction 

Sensitive observations at mid-infrared wavelength regions have been conducted with 
satellites. Infrared Astronomical Satellite (IRAS, Neugebauer et al. 1984), AKARI (Ishihara 
et al. 2010), and Wide-held Infrared Survey Explorer (WISE, Wright et al. 2010) were 
launched in 1983, 2006, and 2009, respectively, and performed all-sky surveys. Infrared Space 
Observatory (Kessler et al. 1996), Midcourse Space Experiment (Price et al. 2001), and Spitzer 
Space Telescope were launched in 1995, 1996, and 2003, respectively, and observed many ob¬ 
jects. Although less sensitive, there are ground-based mid-infrared (MIR) data. 

Using these data, some astronomers have discussed MIR light variations. Meixner et al. 
(2006) conducted Surveying the Agents of a Galaxy’s Evolution (SAGE) of the Large Magellanic 
Cloud with Spitzer Space Telescope. Vijli et al. (2009) detected many near-infrared (NIR) and 
MIR variable objects using the SAGE data taken over 3 months of difference in observing times. 
In the data, a large fraction of the extreme asymptotic giant branch (AGB) stars showed light 
variation. They also detected variable young stellar object (YSO) candidates. Melis et al. 
(2012) reported dramatic MIR flux decrease in young, Sun-like star (TYC 8241 2652 1), using 
AKARI and WISE data and by their own MIR observations with Thermal-Region Camera 
Spectrograph at the Gemini South telescope. They suggested that the reason for the large light 
variation is a disappearance of a warm, dusty circumstellar disk, but how the disk disappeared 
is currently not known. Gandhi, Yamamura, and Takita (2012) combined IRAS, Two Micron 
All-Sky Survey (2MASS) and WISE point source catalogs (PSCs) and found an infrared variable 
source. Based on the photometric characteristics of the object, they concluded that the object is 
likely to be a Sakurai’s object (Duerbeck & Benetti 1996). Rebull (2011) summarized previous 
studies of NIR and MIR variation in young stars and Hoard et al. (2009) reviewed the researches 
of NIR and MIR variation in cataclysmic variables with Spitzer Space Telescope data. 

As described above, some papers that discussed MIR light variations have been pub¬ 
lished. However, the number of papers discussing MIR light variation is much smaller than 
that of papers on optical or NIR light variations because accurate MIR observation is rather 
difficult. We need to observe MIR or far-infrared (FIR) light variation for studying variations in 
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Fig. 1. Normalized response curve. Left: The solid curve corresponds to IRAS 12, the dashed curve to 
AKARI S9W, and the dotted curve to WISE W3. Right: The solid line represents IRAS 25, the dashed 
curve AKARI L18W, and the dotted curve WISE W4. 

circumstellar dust because the peak of dust emission is in MIR or FIR region. Additionally, most 
of the previous works studied short-term light variations (mainly using Spitzer Space Telescope) 
and selected targets by setting some selection criteria. For example, Gandhi, Yamamura, and 
Takita (2012) selected targets by 2MASS and WISE color. 

In this paper, we take advantage of the fact that IRAS, AKARI and WISE performed 
all-sky surveys at similar wavelengths (IRAS 12, AKARI S9W, and WISE W3 and IRAS 25, 
AKARI L18W, and WISE W4, see figure 1). We combined these PSCs and searched the 
objects whose AKARI or WISE MIR fluxes had significantly increased from that of IRAS over 
the 20-30 years of difference in observing times. The aim of this study is to search sources that 
show significant MIR brightenings, and to know what types of objects show MIR large light 
variations and also to know why they show such large light variations. 


2. Data 

2.1. Target Selection 

At first, we merged AKARI/IRC All-Sky Survey Point Source Catalog (Version 1.0, 
Ishihara et al. 2010) and WISE All-Sky Data Release (Cutri et ah 2012) with a tolerance 
angular radius of 5 arcsec. If more than one WISE source were found, we regarded the closest 
one as the identical source. We call it AKARI-WISE catalog. Combining AKARI and WISE 
PSCs is necessary to ensure that the sources are securely detected by both AKARI and WISE 
satellites. Then, we cross-correlated IRAS PSC, Version 2.0 (Helou & Walker 1988) with the 
AKARI-WISE catalog with a tolerance angular radius of 5 arcsec. The beam size of IRAS 
is about 10 times larger than those of AKARI and WISE, making the choice of the tolerance 
radius difficult. If we apply a small radius, the number of the successful cross identification will 
be small. If we use a large radius, we will suffer from miss-identifications. We took 5 arcsec 
as a compromise size between these two problems. When we applied 5 arcsec as a tolerance 
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radius, we did not find the source that had another sources within the tolerance radius. Using 
the IRAS-AKARI-WISE combined catalog, we chose our sample stars by following criteria. 

1. The flux density quality of IRAS at 12 /mi or 25 /irn should be 3 

2. The signal-to-noise ratios (S/N) of IRAS, AKARI and WISE data should be better than 
3 

3. Declination of the sources are larger than —30° so that they can be observed from northern 
hemisphere 


The reason we set the criteria of the flux density quality of IRAS and the S/N of IRAS, AKARI 
and WISE data is to discuss sources with significant flux increasing. The number of sources at 
each step of choosing our sample is shown in table 1. 

The flux density ratios between IRAS and AKARI or WISE versus the numbers of 
sources with these flux density ratio are shown in figure 2, where 

A/ AKARI S9W 


AZ?A12 — 


( 1 ) 
( 2 ) 
( 3 ) 

r u, IRAS 25 
_ A/, WISE W4 

f UtW25 — ~P - (4J 

Uy,IRAS 25 

and F u means flux density. To select only the sources with strongest brightening, we further 
applied the following criteria. 


AAwi2 — 


FRA25 ~ 


Akiras 12 

Fu, WISE W3 

Akiras 12 

F v> AKARI L18W 

Akiras 25 
F u , WISE W4 


1. AR A i2 > 10 

2 . FRyj 12 > 10 

3. FR A25 > VlO ~ 3.16 

4. FR W25 > VlO ~ 3.16 


These numbers were chosen rather arbitrary for this pilot study, so there is no physical meaning 
for the numbers. But as can be seen in figure 2, the above criteria successfully chose sources 
showing the largest brightness changes. As a result, 4 sources were chosen. In the selection 
process, we did not apply color correction. 

In this paper, we do not discuss the sources whose AKARI or WISE’s fluxes decreased 
from that of IRAS, because spatial resolution of IRAS is far worse than that of AKARI or 
WISE. Because of poor spatial resolution, it is possible that IRAS detected some point sources 
as one object which AKARI and WISE were able to resolve. In such cases, IRAS always 
overestimate fluxes compared to AKARI and WISE. 

We chose 4 sources as our targets. The flux density ratios of these sources are shown in 
table 2. Basic information of them is shown in table 3. No other sources are present within 5 arc- 
sec around these 4 sources in IRAS catalog and AKARI-WISE catalog. However, as we show 
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Table 1. The number of sources at each step of choosing our sample. 


Step 

AKARI S9W 

WISE W3 

AKARI L18W 

WISE W4 

1. Merging AKARI and WISE catalogs 

839438 

863972 

191747 

864130 

2. Merging IRAS and AKARI-WISE catalogs 

57973 

61944 

54401 

61938 

3. Choosing sources whose declinations are larger than —30° 

33926 

36545 

31267 

36547 

4. Choosing sources whose S/N are better than 3, 

and IRAS flux density qualities are 3 

32677 

35084 

20588 

22054 


Table 2. Flux density ratios of selected sources. 


Target 

FRA12 

FRw 12 

FRA25 

FRw25 

IRAS 19574+4941 

2.4 ± 0.1 

2.4 ± 0.1 

3.2 ± 0.2 

3.3 ± 0.2 

V2494 Cyg 

25.4 ± 1.6 

54.5 ± 3.3 

7.9 ± 0.4 

11.3 ± 0.6 

IRAS 22343+7501 

3.3 ± 0.6 

12.1 ± 0.5 

1.3 ± 0.2 

3.1 ± 0.2 

V583 Cas 

5.2 ± 0.6 

4.0 ± 0.3 

5.1 ± 0.5 

3.9 ± 0.3 


later, one source (IRAS 22343+7501) is resolved into 4 stars (2MASS J22352345+7517076, 
2MASS J22352442+7517037, [RD95] C, and 2MASS J22352497+7517113) by Cutri et al. 
(2003), Kun et al. (2009), and our observation. These sources are too close to be detected 
individually by IRAS, AKARI, and WISE. We made spectral energy distributions (SEDs) of 
our targets except for IRAS 22343+7501, using all available data to date and also data from 
our observation (see tables 4-9). They are shown in figures 3-5. Meanings of the shapes of 
marks are given in table 10. Errors in observed fluxes are indicated, but it is hard to see these 
errors because they are smaller than the size of marks. Here, we stress that MIR brightening 
is significant. 


Table 3. List of observed stars. 


Target names 

Right Ascension* 

(J2000) 

Declination* 

(J2000) 

IRAS 

AKARI 

WISE 

Other names 

IRAS 19574+4941 

19 58 54.895 

+49 49 49.67 

19574+4941 

1958549+494949 

195854.89+494949.4 


V2494 Cyg 

20 58 21.093 

+52 29 27.70 

20568+5217 

2058211+522927 

205821.09+522927.6 

HH 381 IRS, MSX6C G091.6383+04.4007, WB89 19 

2MASS J22352345+7517076 

22 35 23.458 

+75 17 07.63 

22343+7501 

2235235+751708 

223523.61+751708.3 

SSTc2d J223523.4+751708 

2MASS J22352442+7517037 

22 35 24.425 

+75 17 03.73 

22343+7501 

2235235+751708 

223523.61+751708.3 

e.g. [KP93] 3-8, [RD95] B, SSTc2d J223524.4+751704 

[RD95] C 



22343+7501 

2235235+751708 

223523.61+751708.3 


2MASS J22352497+7517113 

22 35 24.978 

+75 17 11.37 

22343+7501 

2235235+751708 

223523.61+751708.3 

e.g. [RD95] A, [S2006] 8, XMMU J223524.6+751710 

V583 Cas 

23 31 16.203 

+61 32 17.96 

23289+6115 

2331162+613218 

233116.24+613218.3 

CGCS 5887, NIKC 4-136 


Keierence. 


* Cutri et al. (2003) 
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Table 4. Photometric data of IRAS 19574+4941. 


Band 

Wavelength 

A F x 

Epoch 

Reference* 


H 

(Wm -2 ) 



B 

0.44 

1.327 x 10” 14 ± 1.7 x 10 -16 


(1) 

g’ 

0.49 

2 503 x 10~ 14 +2.92x 10 15 
z.ouu x ru — 2.61 xio -15 


(1) 

V 

0.55 

2 215x1 0~ 14 +1-324x10^ 14 
Z.Z1D X 1U — 8.29X10 -15 


(2) 

V 

0.55 

5.186 xio- 14 


(1) 

r’ 

0.63 

9.801 x io- 14 +™:: 6 6 

2001 August 11 

(5) 

r’ 

0.63 

i.026 xio- 13 


(1) 

R 

0.71 

3.723 x IO -13 


(7) 

i’ 

0.77 

6 015 X 10“ 13 +4-84X10- 14 
O.UlO X 1U —4.48xl0~ 15 


(1) 

Ic 

0.81 

1.539 XIO" 12 


(2) 

J 

1.2 

1.069 x IO" 11 

2000 May 30 

(8) 

J 

1.2 

1.025 xlO- 11 + 8 9 2 0x 1 1 0 °: i 1 3 

2013 August 26 

(9) 

H 

1.7 

1 141 x 10~ n + 5 - 0xl0_13 

X 1U -4.8X10- 13 

2000 May 30 

(8) 

H 

1.7 

9.173 xlO" 12 + 9 8 7 4 8 X X 1 1 0 °: 1 1 3 3 

2013 August 26 

(9) 

K s 

2.2 

8.964 xlO" 12 ““ 

2000 May 30 

(8) 

K s 

2.2 

7.938 xlO" 12 + 1 311X 1 1 0 °: 1 1 2 2 

2013 August 26 

(9) 

W1 

3.4 

3 216xl0 -12 + 2 - llxl0 ~ 13 

O.ZIOXIU —1.98xl0~ 13 

2010 

(11) 

W2 

4.6 

2 726 xl0 -12 + 134xl0_13 

Z-'ZDXiU —1.28 X 10 —13 

2010 

(11) 

S9W 

9.0 

2.007 x 10- 12 ±1.1 x 10 -14 

2006-2007 

(15) 

W3 

12 

1.628 x IO” 12 + 8x1 0 °: i I 5 6 

2010 

(11) 

12 

12 

6.442 x 10 _13 ± 2.58 x 10" 14 

1983 

(17) 

L18W 

18 

6.391 x 10 _13 ± 3.4 x IO -15 

2006-2007 

(15) 

W4 

22 

5.373 x IO" 13 + 4 5X 1 1 0 °: 1 1 5 

2010 

(11) 

25 

25 

1.443 x 10 _13 ± 8.7 x 10 -15 

1983 

(17) 

60 

60 

1.998 x 10 -14 (Upper limit) 

1983 

(17) 

100 

100 

4.058 x IO -14 (Upper limit) 

1983 

(17) 


JNotes. 

* Reference numbers are described in table 10. 
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Table 5. Photometric data of V2494 Cyg. 


Band 

Wavelength 

A F x 

Epoch 

Reference* 


(H 

(Wnf 2 ) 



Rc 

0.66 

1.700 x IO" 15 


(6) 

J 

1.2 

9 333 x 10 -14 + 2 - 70xlcr15 
y.ooo x tu _ 2 . 63 xio- 15 

1999 June 21 

(8) 

J 

1.2 

8 979 xl0 -14 + 1114xl0 ~ 14 
s.yiyxiu — 9.91 x io -15 

2013 August 26 

(9) 

H 

1.7 

2 194xl0 -13 + 6 - lxl ° -15 
z.iy^xiu _6.oxio- 15 

1999 June 21 

(8) 

H 

1.7 

9 9QQ -v IO" 13 +2.01xl0 -14 
Z.Z59 X 1U -1.85X10- 14 

2013 August 26 

(9) 

2.05 

2.05 

9.430 x 10" 12 

1996 November 17-19 

(10) 

K s 

2.2 

4 411 x 10 -13 + 7 - 0xl0 ~ 15 

X iU —6.9x 10 -15 

1999 June 21 

(8) 

K s 

2.2 

5 004 x 10 -13 +6-42X10" 14 
0.UU4 X 1U —5.69X10" 14 

2013 August 26 

(9) 

2.50 

2.50 

1.375 x IO" 11 

1996 November 17-19 

(10) 

W1 

3.4 

1 817 x 10“ 12 +i o 7 xio- 13 
1.811 X 1U -l.oixlO- 13 

2010 

(11) 

W2 

4.6 

4 989 x 10 -12 + 3 03xl0 ~ 13 
4.ysy x tu — 2 . 86 xio -13 

2010 

(11) 

IRAC4 

8.0 

3.324 x 10~ 12 ± 3.8 x 10“ 14 

2008 June 29 

(13) 

A 

8.28 

1.382 x 10 -12 ± 5.7 x 10 -14 

1996-1997 

(14) 

S9W 

9.0 

2.908 x 10 -12 ± 3.1 x 10 -14 

2006-2007 

(15) 

W3 

12 

4.854 x 10" 12 ±4.9 x 10~ 14 

2010 

(11) 

12 

12 

8.587 x 10~ 14 ± 6.55 x IO" 15 

1983 

(17) 

C 

12.13 

1.571 x 10 -12 ± 8.3 x 10~ 14 

1996-1997 

(14) 

D 

14.65 

1.750 x 10 _12 ± 1.07 x IO" 13 

1996-1997 

(14) 

L18W 

18 

3.879 x 10 -12 ± 7.3 x 10~ 14 

2006-2007 

(15) 

E 

21.34 

2.550 x 10 _12 ± 1.56 x IO" 13 

1996-1997 

(14) 

W4 

22 

4.506 x 10 -12 ± 2.9 x 10~ 14 

2010 

(11) 

25 

25 

3.533 x 10 _13 ± 1.77 x IO" 14 

1983 

(17) 

60 

60 

4.873 x 10 -13 ± 3.41 x IO" 14 

1983 

(17) 

N60 

65 

3.113 x 10 -12 ± 8.6 x 10 -14 

2006-2007 

(16) 

WIDE-S 

90 

2.023 x 10 _12 ± 1.91 x IO" 13 

2006-2007 

(16) 

100 

100 

3.208 x 10 _13 ± 5.10 x IO" 14 

1983 

(17) 

WIDE-L 

140 

1.273 x 10 _12 ± 1.24 x IO" 13 

2006-2007 

(16) 

N160 

160 

1.116 x 10" 12 ± 6.7 x 10 -14 

2006-2007 

(16) 


Notes. 

* Reference numbers are described in table 10. 
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Table 6. Photometric data of 2MASS J22352345+7517076. 


Band 

Wavelength 

A F x 

Epoch 

Reference* 


(H 

(Wm -2 ) 



J 

1.2 

7.864 x IO -15 (Upper limit) 

1999 October 11 

(8) 

J 

1.2 

3 675 xl0 -15 +8.42xio -16 

O.ulO X 1U —6.85xl0“ 16 

2013 August 26 

(9) 

H 

1.7 

2.848 x 10~ 14 (Upper limit) 

1999 October 11 

(8) 

H 

1.7 

i.oi9 x io- 13 

2013 August 26 

(9) 

K s 

2.2 

2.140 xlO- 14 

1999 October 11 

(8) 

K s 

2.2 

8.281 x IO' 13 

2013 August 26 

(9) 

W1 

3.4 

6 921 x IO -12 + 6 - 54xl0_13 
O.yzi X 1U _5.97xl0- 13 

2010 

(11) 

IRAC1 

3.6 

3.106 x IO -13 ± 2.14 x 10" 14 

2004 April 4 

(12) 

IRAC2 

4.5 

6.662 x IO -13 ± 7.59 x 10" 14 

2004 April 4 

(12) 

W2 

4.6 

2.312 xlO- 11 + 5 5 6 X X 1 1 0 °: 1 1 3 3 

2010 

(11) 

IRAC3 

5.8 

9.614 x IO -13 ± 6.15 x 10" 14 

2004 April 4 

(12) 

S9W 

9.0 

5.495 x 10" 12 ± 1.039 x IO" 12 

2006-2007 

(15) 

W3 

12 

1.556 xlO- 11 + 1 ; 1 2 X X1 0 °: 1 1 3 3 

2010 

(11) 

12 

12 

1.241 x 10 _12 ± 5.0 x 10" 14 

1983 

(17) 

L18W 

18 

5.638 x 10 -12 ± 7.38 x 10" 13 

2006-2007 

(15) 

W4 

22 

1.086 x 10- n ±l x 10" 14 

2010 

(11) 

25 

25 

3.128 x IO -12 ± 1.56 x IO" 13 

1983 

(17) 

60 

60 

3.316 x IO -12 ± 2.98 x IO" 13 

1983 

(17) 

N60 

65 

3.706 x IO -12 ± 4.77 x IO" 13 

2006-2007 

(16) 

MIPS 2 

72 

1.977 x IO -12 ± 4.37 x IO" 13 


(12) 

WIDE-S 

90 

2.415 x IO -12 ± 2.72 x IO" 13 

2006-2007 

(16) 

100 

100 

2.400 x IO -12 ± 3.84 x IO" 13 

1983 

(17) 

WIDE-L 

140 

2.069 x IO -12 ± 1.74 x IO" 13 

2006-2007 

(16) 

N160 

160 

2.024 x IO -12 ± 1.61 x IO" 13 

2006-2007 

(16) 

450 

450 

7.555 x IO -14 


(18) 

800 

800 

2.661 x IO -15 ± 4.27 x IO" 16 

1992 

(19) 

850 

850 

2.003 x IO -14 


(18) 

1100 

1100 

1.044 x 10“ 15 ± 8.2 x IO" 17 

1992 

(19) 

1300 

1300 

5.350 x IO -16 ± 5.53 x IO" 17 

1992 

(19) 


Notes. 

* Reference numbers are described in tabic 10. 


t 2MASS J22352345+7517076, 2MASS J22352442+7517037, [RD95] C, and 2MASS J22352497+7517113 
were not separated without Kun et al. (2009), 2MASS, ISLE, and IRAC (C2D Fall’07 High Reliability 
(HREL) CORES Catalog). 





Table 7. Photometric data of 2MASS J22352442+7517037. 


Band 

Wavelength 

A F x 

Epoch 

Reference* 


W 

(Wm -2 ) 



V 

0.55 

1.281 x 10" 17 

2004-2006 

(3) 

Rc 

0.66 

1.155 xlO- 16 + 2 3 9 7 9 x x 1 1 0 °: i 1 7 7 

2004-2006 

(3) 

Ic 

0.81 

7.140 x 10“ 16 

2004-2006 

(3) 

J 

1.2 

2.540 x 10 -14 (Upper limit) 

1999 October 11 

(8) 

J 

1.2 

1 389 x 10" 14 + 14 4xio _15 

l.ooy X 1U _i.30xl0- 15 

2013 August 26 

(9) 

H 

1.7 

6.744 x 10 -14 (Upper limit) 

1999 October 11 

(8) 

H 

1.7 

4.894 xlO’ 14 + 6 8 5 0 8X 1 1 0 °: 1 1 6 5 

2013 August 26 

(9) 

K s 

2.2 

1.026 xlO’ 13 + 3 3 5 X X 1 1 0 °: 1 1 : 

1999 October 11 

(8) 

K s 

2.2 

8 420xl0 -14 + 9 - 58xl ° 15 

0.4ZU X 1U —S.60x 10 -15 

2013 August 26 

(9) 

W1 

3.4 

6.921 x 10- 12 + 3 ; 34x -::: 

2010 

(11) 

W2 

4.6 

2.312 xlO’ 11 + 5 5 6 X X 1 1 0 °: 1 1 3 3 

2010 

(11) 

S9W 

9.0 

5.495 x 10" 12 ± 1.039 x 10~ 12 

2006-2007 

(15) 

W3 

12 

1.556 xlO’ 11 + 1 ; 1 2X 1 1 0 °: 1 1 3 3 

2010 

(11) 

12 

12 

1.241 x 10 _12 ± 5.0 x 10" 14 

1983 

(17) 

L18W 

18 

5.638 x 10 -12 ± 7.38 x 10" 13 

2006-2007 

(15) 

W4 

22 

1.086 x 10- n ±l x 10" 14 

2010 

(11) 

25 

25 

3.128 x 10 -12 ± 1.56 x 10" 13 

1983 

(17) 

60 

60 

3.316 x 10 -12 ± 2.98 x 10" 13 

1983 

(17) 

N60 

65 

3.706 x 10 -12 ± 4.77 x 10" 13 

2006-2007 

(16) 

MIPS 2 

72 

1.977 x 10 -12 ± 4.37 x 10" 13 


(12) 

WIDE-S 

90 

2.415 x 10 -12 ± 2.72 x 10" 13 

2006-2007 

(16) 

100 

100 

2.400 x 10 -12 ± 3.84 x 10" 13 

1983 

(17) 

WIDE-L 

140 

2.069 x 10 -12 ± 1.74 x 10" 13 

2006-2007 

(16) 

N160 

160 

2.024 x 10 -12 ± 1.61 x 10" 13 

2006-2007 

(16) 

450 

450 

7.555 x lO" 14 


(18) 

800 

800 

2.661 x 10 -15 ± 4.27 x 10" 16 

1992 

(19) 

850 

850 

2.003 x 10" 14 


(18) 

1100 

1100 

1.044 x 10" 15 ±8.2x 10" 17 

1992 

(19) 

1300 

1300 

5.350 x 10 -16 ± 5.53 x 10" 17 

1992 

(19) 


Notes. 

* Reference numbers are described in table 10. 


t 2MASS J22352345+7517076, 2MASS J22352442+7517037, [RD95] C, and 2MASS J22352497+7517113 
were not separated without Kun et al. (2009), 2MASS, ISLE, and IRAC (C2D Fall’07 High Reliability 
(HREL) CORES Catalog). 
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Table 8. Photometric data of 2MASS J22352497+7517113. 


Band 

Wavelength 

A F x 

Epoch 

Reference* 


W 

(Wm -2 ) 



V 

0.55 

8.466 x IO" 17 

2004-2006 

(3) 

Rc 

0.66 

3 180 X 10 -16 +3-224X10- 16 
■J.10U X 1U —1.601 X10 -16 

2004-2006 

(3) 

Ic 

0.81 

2 294 v lfl“ 15 +2-453x10 15 
p.o/n x 1 U — 1.411 x io -16 

2004-2006 

(3) 

J 

1.2 

9.333 xio" 14 

1999 October 11 

(8) 

J 

1.2 

6.547 X 10~ 14 + 6 7 3°6 5 x X iT" 

2013 August 26 

(9) 

H 

1.7 

2 194xl0 -13 + 60xl0 ~ 15 
/.iynxiu _6.0xio- 15 

1999 October 11 

(8) 

H 

1.7 

2 4 8 7 x 10 -13 + 3 - 48xl0 ~ 14 
z.noi x ru _ 3 .o 5 xio- 14 

2013 August 26 

(9) 

K s 

2.2 

4 411xl0 -13 + 70xl0_15 

4.411 X 1U —6.9xl0~ 15 

1999 October 11 

(8) 

K s 

2.2 

2 945 x 10 -13 +4-21 x io~ 14 

,1X140 X 1U -3.80X10- 14 

2013 August 26 

(9) 

W1 

3.4 

6.921 x IO" 12 + 5 6 9 54 X X 1 1 0 0 : 1 1 3 3 

2010 

(11) 

W2 

4.6 

2.312 XIO’ 11 + 5 5 6 x x 1 1 0 °: i 1 3 3 

2010 

(11) 

S9W 

9.0 

5.495 x 10" 12 ± 1.039 x 10~ 12 

2006-2007 

(15) 

W3 

12 

1.556 xlO- 11 + 1 1 ; 1 2X 1 1 0 °: 1 1 3 3 

2010 

(11) 

12 

12 

1.241 x 10 _12 ± 5.0 x 10" 14 

1983 

(17) 

L18W 

18 

5.638 x 10 -12 ± 7.38 x 10" 13 

2006-2007 

(15) 

W4 

22 

1.086 x 10~ 11 ±1 x 10" 14 

2010 

(11) 

25 

25 

3.128 x 10 -12 ± 1.56 x 10" 13 

1983 

(17) 

60 

60 

3.316 x 10 -12 ± 2.98 x 10" 13 

1983 

(17) 

N60 

65 

3.706 x 10 -12 ± 4.77 x 10" 13 

2006-2007 

(16) 

MIPS 2 

72 

1.977 x 10 -12 ± 4.37 x 10" 13 


(12) 

WIDE-S 

90 

2.415 x 10 -12 ± 2.72 x IO" 13 

2006-2007 

(16) 

100 

100 

2.400 x 10 -12 ± 3.84 x IO" 13 

1983 

(17) 

WIDE-L 

140 

2.069 x 10 -12 ± 1.74 x IO" 13 

2006-2007 

(16) 

N160 

160 

2.024 x 10 -12 ± 1.61 x IO" 13 

2006-2007 

(16) 

450 

450 

7.555 x IO" 14 


(18) 

800 

800 

2.661 x 10 -15 ± 4.27 x IO" 16 

1992 

(19) 

850 

850 

2.003 x IO" 14 


(18) 

1100 

1100 

1.044 x 10” 15 ± 8.2 x IO" 17 

1992 

(19) 

1300 

1300 

5.350 x 10 -16 ± 5.53 x IO" 17 

1992 

(19) 


Notes. 

* Reference numbers are described in table 10. 


t 2MASS J22352345+7517076, 2MASS J22352442+7517037, [RD95] C, and 2MASS J22352497+7517113 
were not separated without Kun et al. (2009), 2MASS, ISLE, and IRAC (C2D Fall’07 High Reliability 
(HREL) CORES Catalog). 
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Fig. 2. The histograms of flux density ratios of each band. The empty histograms represent all samples 
selected by 4 steps and the black histograms show our 4 targets (for each step and the number of sources, see 
table 1). The rightmost histograms of each figure include all the sources whose FRai 2 or FR W12 > 10.0, 
or FRa 25 or FRw 25 > 3.2. As can be seen in this figure, objects selected are very rare. 

2.2. Observation 

Little is known about the 4 sources selected in the previous section. For example, we 
looked into the SIMBAD database and no paper was found for 2 sources (IRAS 19574+4941 
and V583 Cas.) To identify what these MIR variable sources are, we carried out JHK S - 
band photometric observations and JHK-band spectroscopic observations from 2013 August 
22 to 28 and on 2014 November 20 with the near-infrared imager and spectrograph ISLE 
(Yanagisawa et al. 2006, 2008), which is attached to Cassegrain focus of the 1.88 m telescope at 
Okayama Astrophysical Observatory. The detector of ISLE is a 1024 x 1024 HgCdTe HAWAII 
array which covers a 4(3 x 4(3 held of view with a pixel scale of 0.254 arcsec/pixel. For 
IRAS 19574+4941, we carried out optical spectroscopic observation at 2014 August 28 with 
Medium And Low-dispersion Long-slit Spectrograph (MALLS), which is attached to Nasmyth 
focus of the NAYUTA 2.0 m telescope at Nishi-Harima Astronomical Observatory. 

We took defocused images for two of our targets (IRAS 19574+4941 and V583 Cas) 
because these stars were too bright. Details of imaging observations for individual stars are 
shown in table 11. 

In observation with ISLE, we obtained both low dispersion (R ~ 510, 350, 410 for J, 
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IRAS 19574+4941 



Fig. 3. The SED of IRAS 19574+4941. Meanings of the shapes of marks are given in table 10. Errors in 
observed fluxes are indicated, but it is hard to see because they are usually smaller than the size of marks. 
Points with arrows represent the data are upper limit. 


V2494Cyg 



Fig. 4. Same as figure 3, but for V2494 Cyg. 
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V583 Cas 



X [pm] 

Fig. 5. Same as figure 3, but for V583 Cas. 

H, K, respectively) and medium dispersion (R ~ 2400, 3600, 2000 for J, H, K, respectively) 
spectra with a slit of 1.0 arcsec (= 4 pixels) width. We also took Ar and Xe comparison 
lamps and dome flats. For medium dispersion mode, we obtained dome flats immediately 
before or after each target exposure to reduce the influence of fringe patterns. We observed A0 
stars to calibrate atmospheric transmissions and instrumental efficiency. In observation with 
MALLS, we obtained low dispersion (R ~ 600) spectra with a slit of 1.2 arcsec width. We also 
took dark frames, Fe-Ne-Ar comparison lamps, dome flats and a spectrophotometric standard 
star. Details of spectroscopic observations of individual stars are shown in table 12. Roughly 
estimated S/N are also listed in the table. 

2.3. Data reduction 

2.3.1. Photometry 

Obtained images were reduced with the standard method, i.e. dark frame subtrac¬ 
tion (including bias-image subtraction), sky subtraction and flat fielding correction. Then 
we combined 8 dithered images by taking average with 3a clipping and performed pho¬ 
tometry for these combined images. We performed aperture photometry for three stars 
(IRAS 19574+4941, V2494 Cyg, and V583 Cas) and point spread function (PSF) fitting pho¬ 
tometry for the others (2MASS J22352345+7517076, 2MASS J22352442+7517037, [RD95] C, 
and 2MASS J22352497+7517113) because they were too crowded to perform aperture photom¬ 
etry. 
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Table 9. Photometric data of V583 Cas. 


Band 

Wavelength 

A F x 

Epoch 

Reference* 


(H 

(Wm -2 ) 



V 

0.55 

8.011 x 10~ 15 


(4) 

Ic 

0.81 

8.274 x 10" 14 


(6) 

J 

1.2 

2 381 x 10~ 13 + 5 - lxl °- 15 

Z.JOl X iU —5.0xl0 -15 

1999 October 15 

(8) 

J 

1.2 

2 158x10~ 13 + 8 - 9xl0 ~ 16 

Z.iOO X iU —8.6xl0 -15 

2013 August 26 

(9) 

H 

1.7 

8.537 x 10- 13 + 1 ; 

1999 October 15 

(8) 

H 

1.7 

8.679 xlO- 13 

2013 August 26 

(9) 

K s 

2.2 

2.138 xlO- 12 

1999 October 15 

(8) 

K s 

2.2 

2 435 x nr 12 +1-77X10- 13 

X 1U -1.65X10- 13 

2013 August 26 

(9) 

W1 

3.4 

o qco v -in —12 +3.20x10 13 
O.OOZ A ±u —2.92x 10 -13 

2010 

(11) 

W2 

4.6 

5 754 x 1 0 -12 +6-3ixio- 13 

°.' 04 X IU —5.69X10” 13 

2010 

(11) 

A 

8.28 

1.152 x 10 _12 ±4.7 x 10 -14 

1996-1997 

(14) 

S9W 

9.0 

1.655 x 10 -12 ± 1.62 x 10" 13 

2006-2007 

(15) 

W3 

12 

9 812xl0 -13 + 3 - 78xl0 ~ 14 

a.oiz x iu -3.64X10- 14 

2010 

(11) 

12 

12 

2.376 x 10 _13 ± 1.43 x 10~ 14 

1983 

(17) 

C 

12.13 

3.981 x 10 _13 ± 2.59 x 10~ 14 

1996-1997 

(14) 

D 

14.65 

2.540 x 10 _13 ± 1.82 x 10" 14 

1996-1997 

(14) 

L18W 

18 

3.016 x 10 _13 ± 1.30 x 10" 14 

2006-2007 

(15) 

W4 

22 

1 894xl0 -13 + 4 - 9xl0_15 

l.oy^ X IU —4.8X10- 15 

2010 

(11) 

25 

25 

4.255 x 10" 14 ±3.40x 10” 15 

1983 

(17) 

60 

60 

1.276 x 10 -13 (Upper limit) 

1983 

(17) 

100 

100 

8.022 x 10~ 13 (Upper limit) 

1983 

(17) 


JNotes. 

* Reference numbers are described in table 10. 


To convert instrumental magnitudes to calibrated apparent magnitudes, we referenced 
the 2MASS All-Sky Catalog of Point Sources (Cutri et al. 2003). We calculated differences 
between 2MASS magnitude and instrumental magnitude of all sources detected in the field of 
view (except the targets) as follows. 

1. We extract instrumental fluxes for a target star and in-held calibration sources (hereafter 
calibrators) in the combined image. 

2. We compute target star’s magnitude for each calibrator using 

o r l „ ( f^X inst,target | /r\ 

^target ~ ™cal = -2.51og 10 — - , (5) 

\ J lux inst, cal / 

where m target is the magnitude of the target star and m ca \ is the magnitude of the chosen 
calibration source from 2MASS. The flux mst, target is the aperture or PSF-fit flux extracted 
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Table 10. Means of point shapes in figures 3, 4, and 5. 


Point shapes 

Reference 

Number 

Plus 

UCAC4 Catalog (Zacharias et al. 2012) 

(1) 

Cross 

TASS Mark IV patches photometric catalog, version 2 (Droege et al. 2006) 

(2) 

Star 

Kun et al. (2009) 

(3) 

Open Box 

Catalog of Stellar Spectral Classifications (Skiff 2013) 

(4) 

Filled box 

Carlsberg Meridian Catalog 14 (CMC14) (Copenhagen University, Obs. et al. 2006) 

(5) 

Open circle 

General Catalog of Variable Stars (Samus et al. 2009) 

(6) 

Filled circle 

Red variables in the NSVS (Wozniak et al. 2004b) 

(7) 

Open triangle 

2MASS All-Sky Catalog of Point Sources (Cutri et al. 2003) 

(8) 

Filled triangle 

ISLE (our observation) 

(9) 

Inverted triangle 

Reipurth & Aspin (1997) 

(10) 

Filled inverted triangle 

WISE All-Sky Data Release (Cutri et al. 2012) 

(11) 

Diamond 

C2D Fall’07 High Reliability (HREL) CORES Catalog 

(12) 

IRS Enhanced Products 

(13) 

Filled diamond 

MSX6C Infrared Point Source Catalog (Egan et al. 2003) 

(14) 

Pentagon 

AKARI/IRC mid-IR all-sky Survey (Ishihara et al. 2010) 

(15) 

AKARI/FIS All-Sky Survey Point Source Catalogues (Yamamura et al. 2010) 

(16) 

Filled pentagon 

IRAS catalogue of Point Sources, Version 2.0 (Helou & Walker 1988) 

(17) 

Barred circle 

Submillimeter-Continuum SCUBA detections (Di Francesco et al. 2008) 

(18) 

quarter-filled circle 

Rosvick et al. (1995) 

(19) 


for the target star while the flux i nst , cal is the same for the chosen calibration source. 
We compute associated uncertainty for each target star magnitude as follows. 


cr, 


m, target ' 


\ 


(Jr 


m, cal 


+ 


2-5 O’flux inst, target | 

,ln 10 flux i ns t , target / 


2-5 a fiux inst:Ca .i 

1 X 


( 6 ) 


yin 10 flux inst , cal / 

where a mt target is the uncertainty on the magnitude of the target star and cr miCa i comes 
from 2MASS for each calibrators. Both afi ux . and (Jfi ux . , consist of two error 

terms. One term is the random error computed by following formula. 

2 5 s — + area x stdev 2 + stdev* 

> qain nsky /* 7 \ 

random, error =-x --—---, (7) 

In 10 flux 

where area means the number of pixels in the aperture, stdev means the standard deviation 
of the sky value, and nsky means the number of pixels used to estimate the sky level (the 
readout noise is negligible). Our targets are very bright, so Poisson noises are dominant in 
random errors and other noises make little contribution. The other term is the error on the 
repeatability of the photometric measurement due to the atmosphere and the instrument. 
We considered the standard deviation of the flux i nst measurements for the target star and 
all calibrators as this error term. Some calibrators are too faint compared to the target star 
and we cannot perform photometry for the calibrators in each individual image. In that 
case, we substituted y/8 x random error for the error, where 8 is the number of images 
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Table 11. Summary of imaging observations. 


Target 

Date 

Band 

Exposure time 

(second x number) 



J 

4x8 

IRAS 19574+4941 

2013 Aug 26 

H 

4x8 



K s 

4x8 



.1 

4x8 

V2494 Cyg 

2013 Aug 26 

H 

4x8 



K s 

4x8 



J 

30 x 8 

2MASS J22352345+7517076 

2013 Aug 26 

H 

4x8 



K s 

4x8 



J 

30 x 8 

2MASS J22352442+7517037 

2013 Aug 26 

H 

4x8 



K s 

4x8 



J 

30 x 8 

[RD95] C 

2013 Aug 26 

H 

4x8 



K s 

4x8 



J 

30 x 8 

2MASS J22352497+7517113 

2013 Aug 26 

H 

4x8 



K s 

4x8 



J 

4x8 

V583 Cas 

2013 Aug 26 

H 

4x8 



K s 

4x8 


combined. We summed the two terms in quadrature. 

4. We average all target star magnitudes, then average all errors. Three iterations of 3cr 
clipping were used to remove outliers. The number of 2MASS sources in each individual 
image and the number of calibrators are listed in table 13. 

2.3.2. Spectroscopy 

We performed the standard data reduction for all obtained spectra, including sky sub¬ 
traction (which contained bias and dark frame subtraction), flat fielding, and wavelength cali¬ 
bration for NIR spectra. Flat fielding also has a role to de-fringe for medium dispersion spectra. 
Although we performed flat-field correction, fringe features are seen in H-band spectrum of V583 
Cas. Then we tried Fourier technique to de-fringe, but unfortunately, it did not work well. To 
correct the atmospheric transmissions and the instrumental efficiency, spectra of our targets 
were divided by spectra of AO-type stars (HR 7734, HR 7782, HR 8246, HR 8598, HR 8844, 
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Table 12. Summary of spectroscopic observations. 


Target 

Date 

Dispersion 

Band 

Exposure time 

Standard star 

minimum S/N* 





(second x number) 




2014 Aug 28 

low 

optical 

300 x 9 

V2011 Cyg 

29 




J 

20 x 4 


69 



low 

HK 

20 x 4 

HR 7782 

112 

IRAS 19574+4941 

2013 Aug 22 


J 

30 x 6 


141 



medium 

H 

30 x 4 

HR 7734 

77 




K 

30 x 4 


132 




H 

60 x 20 


46 

V2494 Cyg 

2013 Aug 26 

medium 

K 

60 x 10 

HR 8246 

118 


2013 Aug 25 

low 

HK 

30 x 10 

HR 8844 

58 

2MASS J22352345+7517076 

medium 

K 

120 x 4 

HR 8598 

161 

2013 Aug 28 


2013 Aug 25 

low 

HK 

120 x 8 


106 

2MASS J22352442+7517037 

medium 

K 

120 x 20 

HR 8844 

118 

2013 Aug 28 


2013 Aug 25 

low 

HK 

60 x 10 

HR 8844 

151 

2MASS J22352497+7517113 

medium 

K 

120 x 6 

HR 8598 

138 

2013 Aug 28 




J 

120 x 4 


29 



low 

HK 

30 x 4 


100 


2013 Aug 22 


J 

120 x 10 


28 

V583 Cas 

medium 

H 

60 x 12 

HR 9019 

42 





K 

30 x 4 


79 


2014 Nov 20 

low 

HK 

60 x 4 


134 


JNotes. 


* Without the regions referred in captions for figures 6-11. 


and HR 9019) that were observed at similar airmass. These AO-type stars were observed imme¬ 
diately before or after each target exposure. When we performed this correction, the hydrogen 
absorption features (Paschen series lines and Brackett series lines) of each star were removed 
by spectral fitting with Gaussian profile. We assumed that all AO-type stars have spectra of 
10000 K black body. For the optical spectrum, we also performed the standard data reduction, 
which includes dark frame subtraction, flat fielding, sky subtraction, and wavelength calibra¬ 
tion. We used a spectrophotometric standard star (V2011 Cyg) to correct the atmospheric 
transmissions and the instrumental efficiency. 

The dispersion of obtained spectra is 5.9 A/pixel, 13.4 A/pixel, 13.4 A/pixel for J, 
H, K, respectively (low dispersion mode), 1.66 A/pixel, 1.63 A/pixel, 3.36 A/pixel for J, H, 
K, respectively (medium dispersion mode), and 2.89 A/pixel for optical data. The accuracy 
of wavelength determination in wavelength calibration is better than 1.3 A in low dispersion 
mode of ISLE, 0.18 A in medium dispersion mode of ISLE, and 0.2 A for MALLS’s data. 
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Table 13. The number of 2MASS sources seen in each image and 2MASS sources used for photometric calibration. 


Target 

Band 

All 2MASS sources 

2MASS sources used for photometric calibration 


J 

7 

7 

IRAS 19574+4941 

H 

3 

3 


K s 

3 

2 


J 

7 

7 

V2494 Cyg 

H 

9 

9 


K s 

10 

10 


J 

12 

10 

2MASS J22352345+7517076 

H 

13 

10 


K s 

14 

10 


J 

12 

10 

2MASS J22352442+7517037 

H 

13 

10 


K s 

14 

10 


J 

12 

10 

[RD95] C 

H 

13 

10 


K s 

14 

10 


J 

12 

10 

2MASS J22352497+7517113 

H 

13 

10 


K s 

14 

10 


J 

8 

8 

V583 Cas 

H 

11 

11 


K s 

8 

8 


3. Results and Discussions 

The results of photometric calibration are listed in table 14. The reduced spectra of our 
targets are shown in figures 6-12 and features seen in the spectra are listed in tables 15 and 16. 
We discuss each star below. 

3.1. IRAS 19574+4941 

This star showed significant but relatively weak MIR light variations among our targets. 
This star is classified into a semiregular or an irregular variable by Wozniak et al. (2004b) 
using NS VS (Wozniak et al. 2004a) data, and an irregular variable by AAVSO International 
Variable Star Index VSX (Watson et al. 2013). The light curve of this star from NSVS catalog 
is shown in figure 13. 

This star showed no discernible NIR light variations from 2MASS data. This star showed 
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1.10 1.15 1.20 1.25 1.30 

A, [|im] 

Fig. 6. Low dispersion spectra of our targets in the J-band. For comparison purpose, M-type and carbon 
stars are also shown from the IRTF Spectral Library (Rayner et al. 2009). The region of strong atmospheric 
absorption is shown in dark gray. The light gray hatch represents the region where we may not sufficiently 
correct hydrogen absorption features in the calibration process. Stars with asterisks are our target stars. 
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X [|iin] 


Fig. 7. Low dispersion spectra of our targets in the H-band. For comparison purpose, M-type and carbon 
stars are also shown from the IRTF Spectral Library (Rayner et al. 2009). The light gray hatch represents 
the region where we may not sufficiently correct hydrogen absorption features in the calibration process. 
Stars with asterisks are our target stars. 
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Fig. 8. Low dispersion spectra of our targets in the K-band. For comparison purpose, M-type and carbon 
stars are also shown from the IRTF Spectral Library (Rayner et al. 2009). The light gray hatch represents 
the region where we may not sufficiently correct hydrogen absorption features in the calibration process. 
Stars with asterisks are our target stars. 
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Fig. 9. Medium dispersion spectra of our targets in the J-band. For comparison purpose, M-type and 
carbon stars are also shown from the IRTF Spectral Library (Rayner et al. 2009). The region of strong 
atmospheric absorption is shown in dark gray. Stars with asterisks are our target stars. 
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A [|iin] 

Fig. 10. Medium dispersion spectra of our targets in the H-band. For comparison purpose, M-type and 
carbon stars are also shown from the IRTF Spectral Library (Rayner et al. 2009). The light gray hatch 
represents the region where we may not sufficiently correct hydrogen absorption features in the calibration 
process. Features seen in V583 Cas spectrum longer than 1.7 /iin may be fringe features that we cannot 
sufficiently correct. Stars with asterisks are our target stars. 
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Fig. 11. Medium dispersion spectra of our targets in the K-band. For comparison purpose, M-type and 
carbon stars are also shown from the IRTF Spectral Library (Rayner et al. 2009). The light gray hatch 
represents the region where we may not sufficiently correct hydrogen absorption features in the calibration 
process. Stars with asterisks are our target stars. 
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jOw dispersion optical spectra of our targets. For comparison purpose, M-type stars are shown 
.es (1998) and solar metallicity model spectra of carbon stars are shown from Aringer et al. 
Star with an asterisk is our target star. 
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Table 14. Photometric data of our targets. 


Target 

Band 

ISLE mag 

ISLE mag uncertainties 

2MASS mag 

2MASS uncertainties 


J 

6.442 

0.091 

6.397 

0.026 

IRAS 19574+4941 

H 

5.760 

0.110 

5.523 

0.047 


K s 

5.167 

0.166 

5.035 

0.026 


J 

11.586 

0.127 

11.544 

0.031 

V2494 Cyg 

H 

9.770 

0.092 

9.813 

0.030 


K s 

8.168 

0.131 

8.305 

0.017 


J 

15.056 

0.224 

14.230 

Upper limit 

2MASS J22352345+7517076 

H 

10.646 

0.138 

12.030 

Upper limit 


K s 

7.621 

0.126 

11.590 

0.054 


J 

13.612 

0.107 

12.957 

Upper limit 

2MASS J22352442+7517037 

H 

11.442 

0.137 

11.049 

Upper limit 


K s 

10.103 

0.117 

9.888 

0.036 


J 

14.365 

0.106 



[RD95] C 

H 

12.781 

0.132 




K s 

11.889 

0.197 




J 

11.929 

0.111 

11.928 

0.031 

2MASS J22352497+7517113 

H 

9.677 

0.142 

9.714 

0.031 


K s 

8.426 

0.110 

8.461 

0.028 


J 

10.634 

0.044 

10.527 

0.023 

V583 Cas 

H 

8.320 

0.051 

8.338 

0.023 


K s 

6.450 

0.076 

6.591 

0.018 


TiO absorption features in optical and J-band spectra and CO absorption features in H- and 
K-band. These features are typical of an M-type star’s spectrum. Additionally, we calculated 
reduced chi-squares by comparing this star’s spectra with comparison stars’ spectra or model 
spectra (see table 17). When we calculated the reduced chi-squares, we interpolated the values 
of flux densities using cubic spline. The reduced chi-squares were calculated by the following 
formula, 


X 2 = 


N- 1 


N 

E 


(XF x , 


target, i 


XF. 


i” A, comparison or model 


,0 s 


XF, 


( 8 ) 

i = \ /x i M A, comparison or model, i 

where N is the number of data points, Fy target is the flux density of a target’s spectrum, and 
Fa, comparison or model is the flux density of a comparison star’s or model spectrum. The same 
calculation is done for V583 Cas’s spectra. These results suggest that this star is a relatively 
early M-type star. Therefore, we conclude that this star is an M-type star. 

Can this star be a nearby M-dwarf? According to Cox (2000), absolute magnitude of 
the MOV is 5.99, 5.32, and 5.15 for J-, H-, and K-band respectively. By comparing these 
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Table 15. Spectral features seen in each star. 


Target 

Wavelength (/im) 

Features 


0.622 

TiO 


0.705 

TiO 


0.767 

TiO 


0.844 

TiO 


1.104 

TiO 


1.242 

TiO 


1.558 

12 CO 


1.578 

12 CO 


1.598 

12 CO 


1.619 

12 CO 


1.640 

12 CO 


1.662 

12 CO 

IRAS 19574+4941 

1.684 

12 co 


1.707 

12 co 


1.711 

Mg i 


1.723 

12 co 


1.754 

12 co 


2.294 

12 co 


2.323 

12 co 


2.335 

Na i 


2.339 

Na i 


2.345 

13 co 


2.353 

12 co 


2.374 

13 co 


2.383 

IsO 

o 

o 


1.662 

12 co 


1.684 

12 co 


1.707 

12 co 


1.711 

Mg 1 

V2494 Cyg 

2.294 

12 co 

2.323 

12 co 


2.345 

13 co 


2.353 

12 co 


2.374 

13 co 


2.383 

12 co 
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Table 16. Continuation of table 15. 


Target 

Wavelength (//m) 

Features 


2.294 

12 CO 


2.323 

12 CO 


2.345 

13 CO 

2MASS J22352442+7517037 

2.353 

12 CO 


2.374 

13 CO 


2.383 

12 CO 


2.294 

12 CO 


2.323 

12 CO 


2.345 

13 CO 

2MASS J22352497+7517113 

2.353 

12 CO 


2.374 

13 CO 


2.383 

12 CO 


1.53 

HCN + C 2 H 2 


1.77 

c 2 


2.294 

12 co 

V583 Cas 

2.323 

12 co 


2.353 

12 co 


2.383 

12 co 


values with the 2MASS magnitude (table 14), we can estimate its distance from the Earth, 
and it is ~10 pc. However, the proper motion of this star is very small (~10 mas/yr, Monet 
et al. 2003; Roeser et al. 2010), and this contradicts any close distance. Therefore, the star 
is likely to be more distant and hence could be an evolved star. Additionally, this star is not 
a member of known star forming regions. Therefore, this star is probably an AGB star. Few 
monitoring observations for AGB stars in MIR region have been conducted to date. From these 
observations, the amplitudes of Mira variables in MIR are less than 1 mag (Smith et al. 2002). 
However, this star shows stronger brightening than 1 mag at 25 /iin region, so we cannot 
explain the reason for brightening just by pulsation. We infer that a large amount of mass loss 
of the central star is one possible reason for the brightening in MIR. If mass loss rate increases 
an order of magnitude, we can explain brightening to several times, as we assume symmetric 
mass loss. The other possibility is that this star can be in a binary system. The light curve of 
this star in optical is similar to V Hya (Knapp et al. 1999) in terms of having two variability 
periods. V Hya is considered to be in a binary system. Knapp et al. (1999) suggested that 
the shorter period of V Hya (~530 days) was due to pulsation and the longer period (~6000 
days) was due to extinction by circumstellar dust orbiting with the companion. Although the 








Fig. 13. The light curve of IRAS 19574+4941. 

periods of these two stars are different (the periods of IRAS 19574+4941 are ~60 days and 500 
days), we can consider the possibility that this star is also in a binary system. Also, if this 
star is in a binary system, we can consider other reasons for brightening in MIR. One possible 
thing is eclipse, but it has a problem. If the longer period is due to eclipse, the amplitude 
in R-band (~0.7 mag) is smaller than the brightening in MIR. Therefore, eclipse would not 
become main factor of brightening. Another possible thing is mass transfer. When gas and 
dust accrete to companion, gravitational energy is released in the form of light. As a result, 
dust around the system absorbs the light and increases the amount of reradiation. Anyway, 
it needs to confirm whether this star is in a binary system or not to explore the possibility of 
these things. However, it is difficult to confirm whether a red giant is in a binary system or 
not because we cannot always distinguish orbital motion and pulsation. If the system has an 
accretion disk, then ultraviolet spectroscopy should reveal excess continuum or line emission 
from the accretion shock. 

3.2. V2494 Cyg 

This star showed the strongest brightening in MIR in our samples. Interestingly, 
strong brightening was also seen in FIR (F u akari mo/F Ut iras 60 = 6.9 ± 0.5 and 
F v . AKARI WIDE-s/ Fv, IRAS 100 = 5.7 ± 0.8). 

From previous works, this star is considered a FU Orionis star (FUor, e.g. Reipurth 
& Aspin 1997; Greene et al. 2008). A FUor is a kind of YSO, which exhibit violent outburst 
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Table IT. Reduced chi-square values calculated for our targets and comparison stars. 


Star 

Dispersion 

Band 

IRAS 19574+4941 

V583 Cas (2013/8/22) 

V583 Cas (2014/11/20) 



optical 

68.63 





J 

3.922 

27.60 



low 

H 

3.333 

21.10 

29.00 

MOIII 


K 

1.272 

40.63 

55.51 



J 

0.783 

2.340 



medium 

H 

1.860 





IC 

5.181 

16.73 




optical 

33.08 





J 

3.783 

25.24 



low 

H 

2.126 

18.78 

25.96 

M4III 


IC 

1.667 

49.27 

67.11 



.1 

0.655 

2.166 



medium 

H 

1.956 





IC 

4.873 

21.18 




optical 

56.23 





J 

42.95 

12.00 



low 

H 

102.5 

9.149 

9.871 

M8III 


K 

23.29 

11.50 

14.47 



J 

13.87 

2.243 



medium 

H 

17.68 





K 

15.04 

3.229 


C5 (4000 K) 

low 

optical 

82.91 



C7 (3500 IC) 

low 

optical 

77.31 



C9 (3000 IC) 

low 

optical 

73.81 





.1 

9.050 

10.62 



low 

H 

26.80 

10.27 

12.11 



IC 

3.817 

44.06 

59.07 

C-N4.5 


.1 

4.234 

2.828 



medium 

H 

6.721 





IC 

4.701 

18.76 




J 

8.348 

10.28 



low 

H 

39.16 

4.626 

5.197 

C7,6e(N4) 


IC 

6.747 

24.67 

32.87 


J 

2.156 

1.695 



medium 

H 

3.267 





K 

7.755 

16.15 
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because of a sudden increase of the mass accretion rate (up to 10“ 4 M 0 ?/r _1 ) from the accretion 
disk around it (Hartmann & Kenyon 1985). During the outbursts, these stars can increase the 
bolometric luminosity by two to three orders of magnitude. 

This star showed no discernible NIR light variations from 2MASS data. No variation 
in spectrum of this star from previous data (e.g. Reipurth & Aspin 1997; Greene et al. 2008) 
was seen. The deep first overtone CO absorption bands at 2.29 //m were seen, which is usually 
observed in the spectrum of FUor type star. 

Several photometric monitorings for YSOs in MIR have been conducted to date. 
Morales-Calderon et al. (2011) monitored YSOs in Orion Nebula Cluster and Flaherty et al. 
(2012) surveyed five evolved protoplanetary disks in the IC 348 cluster. They considered the 
reason for MIR variability as mass accretion, flares, photospheric spots, or variable extinction. 
The amplitudes of these light variations are smaller than 1 mag and much smaller than that 
of V2494 Cyg. Muzerolle et al. (2013) reported multi-epoch MIR observations of protostar 
LRLL 54361 that exhibits by far the largest MIR flux variability (~2.5 mag) in IC 348. This 
protostar showed periodic MIR flux variation caused by pulsed accretion. Large light variation 
in MIR is also reported by Fischer et al. (2012). They reported multiwavelength observations 
of V2775 Ori, a FUor type star, and strong brightening was also caused by increasing mass 
accretion. To cause such large light variation, it needs to change accretion rate dramatically. 
Therefore, we think the reason MIR flux was significantly increased is a sudden increase of the 
mass accretion rate. It caused brightening of the central star and dust around the star was 
heated rapidly. Heated dust emits more MIR flux. Magakian et al. (2013) suggested that the 
increase in brightness was started in the early 1980s. Our result indicates that this star started 
strong brightening between the time when IRAS and AKARI data were taken, restricting the 
epoch of brightening stronger (after 1983). 

Significant light variation of FUors in optical wavelength was well known, but large light 
variation of MIR and FIR is rarely known. FUors are rarely detected because the brightening 
timescale of FUors (from a few years to several decades) is much shorter than the timescale 
between outbursts of the order of 10 4 yr. 

3.3. 2MASS J2235234 5+7517076, 2MASS J22352U2+ 7517037, [RD95] C, and 

2MASS J223524-97+7517113 

WISE W3 data showed strong brightening from IRAS 12 data. 

These stars are in Cepheus flare region, a star-forming region. Kun et al. 
(2009) carried out optical spectroscopic observations for 2MASS J22352442+7517037 and 
2MASS J22352497+7517113. They classified 2MASS J22352442+7517037 into an MO-type 
star, 2MASS J22352497+7517113 into a K7-type star, and both stars into classical T Tauri 
stars. 

2MASS J22352345+7517076 showed significant light variations in H and K s -band. 
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Especially, it changed about 4 magnitudes in K s -band but usual AGB stars do not show such 
large light variation in NIR (Ita et ah, in prep.). It supports the idea that this star is a 
YSO. The spectrum of this star is nearly featureless and we can barely see CO absorption 
features in K-band. It is probably due to significant excess emission from a disk of material 
that veils the stellar features. It also supports that this star is young star with a protoplane¬ 
tary disk. Therefore, we infer that this star is probably a YSO. 2MASS J22352442+7517037 
and 2MASS J22352497+7517113 have CO absorption features. It is consistent with Kun et al. 
(2009). 

Dramatic brightening is seen in only WISE data, so this brightening occurred between 
2006 and 2009. Because 2MASS J22352345+7517076 showed significant light variations, this 
star probably related to the brightening in MIR. However, IRAS, AKARI, and WISE were 
not able to detect 2MASS J22352345+7517076, 2MASS J22352442+7517037, [RD95] C, and 
2MASS J22352497+7517113 separately. For more detailed discussion, we need the monitoring 
data that can resolve these stars in MIR. 

3.4. V583 Cas 

This star became ~ 4-5 times brighter in MIR. 

This star is classified into LB-type by AAVSO International Variable Star Index VSX 
and General Catalogue of Variable Stars (GCVS, Samus et al. 2009). The spectral type of this 
star is classified into M-type by Skiff (1997), GCVS, and Skiff (2013). A search through the 
literature showed that their spectral classification is based on Rosino, Bianchini, & di Martino 
(1976). On the other hand, Ita et al. (2010) classified this star into carbon star based on its 
red color. Both Kurtanidze & Nikolashvili (1989) and Alksnis et al. (2001) list this star as a 
carbon star. However, their data are not publicly available, so it is not possible to confirm the 
spectral type from previous classifications. 

This star also showed no discernible light variations from 2MASS data. In the H-band 
spectrum of this star, we can see C 2 Ballik-Ramsay system around 1.77 /irn. It is typical of a 
carbon star’s spectrum. In addition to this feature, we can see a weak feature around 1.52 /irn. 
According to Joyce (1998), some carbon stars have an absorption band at 1.53 /mi due to HCN 
and C 2 H 2 . Spectra of these stars show greatly weakened CO absorption features and these 
stars are very red. These characteristics are consistent with this star ((J — iV,) 2 mass = 4.215, 
(J —JC)isle = 3.936). It also suggests that this star is a carbon star. Hence, we conclude that 
this star is a carbon star. 

We cannot see C 2 Ballik-Ramsay system in the spectrum taken on 2013 August 22. The 
reason is probably that absorption by terrestrial water vapor is stronger in summer than in 
autumn in Japan. The wavelength of the feature is near band end of the H-band filter, so S/N 
tends to be worse compared to the other wavelengths. Therefore, C 2 Ballik-Ramsay system 
could have been obscured by noise, and that is why we could not see the feature. 
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As described for IRAS 19574+4941, we cannot explain the significant brightening just 
by pulsation. Therefore, additional mechanisms are needed to explain the MIR brightening. 
However, observational data for this star are insufficient and we cannot confirm the reason for 
brightening now. 

4. Summary 

Using IRAS, AKARI, and WISE point source catalogs, we found that 4 sources 
(IRAS 19574+491, V2494 Cyg, IRAS 22343+7501, and V583 Cas) significantly brightened 
at MIR wavelengths over the 20-30 years of difference in observing times. IRAS 22343+7501 
consists of 4 stars (2MASS J22352345+7517076, 2MASS J22352442+7517037, [RD95] C, and 
2MASS J22352497+7517113) and we obtained first JHK S photometric data for all 4 sources. 
Our spectroscopic observation reveals that two stars among the four are evolved stars (one M- 
type star and one carbon star). The other sources except [RD95] C (not known the type yet) are 
young stars (one FUor, two classical T Tauri star, and one possible YSO). IRAS 19574+4941 
and V583 Cas are especially interesting stars because we cannot explain the reason for signif¬ 
icant brightening in MIR just by pulsation. We infer that large amounts of mass loss of the 
central star and/or phenomena related to a binary star system are possible reasons for the 
brightening in MIR. We need further observations to confirm these inferences. 
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